Controlling Spin Exchange Interactions of Ultracold Atoms in Optical Lattices 
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We describe a general technique that allows to induce and control strong interaction between spin 
states of neighboring atoms in an optical lattice. We show that the properties of spin exchange 
interactions, such as magnitude, sign, and anisotropy can be designed by adjusting the optical 
potentials. We illustrate how this technique can be used to efficiently "engineer" quantum spin 
systems with desired properties, for specific examples ranging from scalable quantum computation 
to probing a model with complex topological order that supports exotic anyonic excitations. 



PACS numbers:03.75.Fi,03.67.-a,71.10.-w,42.50.-p 
Recent observations of the superfluid to Mott insulator 
transition in a system of ultracold atoms in an optical lat- 
tice open fascinating prospects for studying many body 
phenomena associated with strongly correlated systems 
in a highly controllable environment [1-4]. For instance, 
the recent studies have shown that, with spinor bosonic 
or fermionic atoms in optical lattices, it may be possible 
to observe complex quantum phase transitions [5,6], to 
probe novel superfluidity mechanisms [7,8], or to demon- 
strate the spin-charge separation predicted from the Lut- 
tinger liquid model [9]. 

This Letter describes a general technique to control 
many-body spin Hamiltonians using ultra-cold atoms. 
Specifically, we show that when two-state bosonic or 
fermionic atoms are confined in an optical lattice, the in- 
teraction between spins of the particles can be controlled 
by adjusting the intensity, frequency, and polarization of 
the trapping light. The essential idea is to induce and 
control virtual spin-dependent tunneling between neigh- 
boring atoms in the lattice that results in a controllable 
Heisenberg exchange interaction. By combining this sim- 
ple experimental technique with the design of the lattice 
geometry, it is possible to "engineer" many interesting 
spin Hamiltonians corresponding to strongly correlated 
systems. 

Such techniques are of particular significance since 
quantum magnetic interactions are central to under- 
standing complex orders and correlations [10]. We il- 
lustrate this with several examples: (i) we show that one 
of the generated Hamiltonians provides us an easy way 
to realize the so-called cluster states in two or three di- 
mensions [11], which are useful for an implementation 
of scalable quantum computation with neutral atoms; 
(ii) we show that the realized Hamiltonian has a rich 
phase diagram, opening up the possibility to observe vari- 
ous quantum magnetic phase transitions in a controllable 
way; (iii) finally, we show how to implement an exactly- 
solvable spin Hamiltonian recently proposed by Kitaev 
[12], which supports abelian and non-abelian anyonic ex- 
citations with exotic fractional statistics. Abelian anyons 
could also exist in a fast rotating condensate [13]. 

We consider an ensemble of ultracold bosonic or 
fermionic atoms confined in an optical lattice formed by 
several standing wave laser beams. We are interested 



in the Mott insulator regime, and the atomic density of 
roughly one atom per lattice site. Each atom is assumed 
to have two relevant internal states, which are denoted 
with the effective spin index a =t,|, respectively We 
assume that the atoms with spins a =|, J, are trapped by 
independent standing wave laser beams through polar- 
ization (or frequency) selection. Each laser beam creates 
a periodic potential V llr7 sin 2 (k fl ■ r) in a certain direc- 
tion n, where fc^ is the wavevector of light. For suffi- 
ciently strong periodic potential and low temperatures, 
the atoms will be confined to the lowest Bloch band as 
has been confirmed from experiments [1], and the low 
energy Hamiltonian is then given by 
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Here 



denotes the near neighbor sites in the direction 
fx, are bosonic (or fermionic) annihilation operators 
respectively for bosonic (or fermionic) atoms of spin a 
localized on site i, and 7ij CT = a\ a^. 

For the cubic lattice (/i = x,y,z) and using a 
harmonic approximation around the minima of the 
potential [3], the spin-dependent tunneling energies 
and the on-site interaction energies are given by 
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Here V^ n = 
is the spin average poten- 
Er = h 2 k 2 /2m is the atomic 
a s n is the scattering length be- 
spins. For bosonic 
(a sa are 



tial in each direction 
recoil energy, and 
tween the atoms of different 
atoms U a ~ (8/ir) 1/2 (ka sa ) (E R V la V 2 aV 3a 
the corresponding scattering lengths). For fermionic 
atoms, U a is on the order of Bloch band separa- 
tion ~ 2yj / V fJ-a E R , which is typically much larger 
than fjfj and can be taken to be infinite. In writ- 
ing Eq. (1), we have neglected overall energy shifts 

{V E R V rt ~ V E R V vl) Kt ~ n a) A which can be 
easily compensated by a homogeneous external magnetic 

field applied in the z direction. 

From the above expressions, we observe that t^a de- 
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pend sensitively (exponentially) upon the ratios V^/Er 
while Ufi and U a exhibit only weak dependence. So 
we can easily introduce spin-dependent tunneling t^ a by 
varying the potential depth and V^i with control of 
the intensity of the trapping laser. We now show that this 
simple experimental method provides us a powerful tool 
to engineer many-body Hamiltonians. We are interested 
in the regime where <C U a , Ufj and (n^) + (n^) ~ 1, 
which corresponds to an insulating phase. In this regime, 
the terms proportional to tunneling t^ a can be consid- 
ered via perturbation theory. We use a simple generaliza- 
tion of the Schriffer-Wolf transformation [14] (see another 
method in [8]), and to the leading order in t^/U^i, Eq. 
(1) is equivalent to the following effective Hamiltonian 

# =E ± 'V-L (afaf + afa])} , (2) 

(id) 

Here of = - n^, of = a^a^ + a}^a^, and of = 
—i ^a^fti| — a\^a,i^j arc the usual spin operators. The 

+ and — signs before X^± in Eq. (4) correspond respec- 
tively to the cases of fcrmionic and bosonic atoms. The 
parameters A MZ and A M j_ for the bosonic atoms are given 
by 



A 



j-2 — I — /2 



t 



fiZ 



2U 



11 



17 T 



Ui 



11 



(3) 



For fermionic atoms the expression for A^ is the same as 
in (3), but in the expression for A z the last two terms 
vanish since U a U^. In writing Eq. (2), we neglected 

the term 4 (t 2 ^/U^ — t^JU^ erf, which can be eas- 
ily compensated by an applied external magnetic field. 

The Hamiltonian (2) represents the well-known 
anisotropic Heisenberg model (XXZ model), which 
arises in the context of various condensed matter systems 
[10]. However, the approach involving ultracold atoms 
has a unique advantage in that the parameters X^ z and 
\fj_± can be easily controlled by adjusting the intensity 
of the trapping laser beams. They can also be changed 
within a broad range by tuning the ratio between the 
scattering lengths a s ^i and a SIJ (er =f , J,) by adjusting an 
external magnetic field through Feshbach resonance [15]. 
Therefore, even with bosonic atoms alone, it is possible 
to realize the entire class of Hamiltonians in the general 
form (2) with an arbitrary ratio A^z/A^j,. This is impor- 
tant since bosonic atoms are generally easier to cool. In 
Figure la we show the phase diagram of the Hamiltonian 
(2) on a bipartite lattice as a function of fit = Al+^l Al 
and Un/U a [16] for the case when U-\ = U~i and t^ a in- 
dependent of the spatial direction [i. Certain lines on 
this phase diagram correspond to well known spin sys- 
tems: when Ufi/U a = 1/2 we have an XY model; when 
fit = oo (if or t± is zero) we have an Ising model; for 
fit = ±(1/2 - U^i/U^y 1 we have an SU(2) symmetric 
antiferromagnetic or ferromagnetic systems, respectively. 



Before proceeding, we estimate the typical energy 
scales and discuss the influence of imperfections and 
noise. For Rb atoms with a lattice constant it/ k ~ 
426nm, the typical tunnelling rate t/h can be chosen from 
zero to a few kHz [1]. The on-site interaction U/% corre- 
sponds to a few kHz at zero magnetic field, but can be 
much larger near the Feshbach resonance. The energy 
scale for magnetic interaction is about t 2 /TiU ~ 0.1kHz 
(corresponding to a time scale of 10ms) with a conser- 
vative choice of U ~ 2kHz and (t/U) 2 ~ 1/20. These 
energy scales are clearly compatible with current exper- 
iments [1]. We further note that the present system 
should be quite robust to realistic noise and imperfec- 
tions. First of all, the next order correction to the Hamil- 
tonian (2) is proportional to (t/U) 2 , which is small in 
the Mott regime. Second, since the atoms only virtu- 
ally tunnel to the neighboring sites with a small prob- 
ability (t/U) 2 , the dephasing rate and the inelastic de- 
cay rate are significantly reduced compared with the cold 
collision scheme [17,20]. Finally, the spontaneous emis- 
sion noise rate can be made very small by using a blue- 
detuned optical lattice or by increasing the detuning. In 
a blue-detuned lattice, even with a moderate detuning 
A ~ 5GHz, the effective spontaneous emission rate is es- 
timated to be of the order of Hz, which is significantly 
smaller than t 2 / (TiU). 

We now illustrate the ability to engineer many-body 
spin Hamiltonians with specific examples. For the first 
example, we set V^/V^ 3> 1, so that t^i becomes negli- 
gible while t M | remains finite. In this case, the Hamilto- 
nian (2) reduces to the Ising model H = J2(ij) ^t^z^i^j, 
with A M 2 = t 2 ^/ (0.5/U n - 1/Uf). Though this Hamilto- 
nian has quite trivial properties for its ground states and 
excitations, its realization in optical lattices can be very 
useful for a dynamical generation of the so-called cluster 
states [11]. Specifically, we note that this Ising interac- 
tion can be easily turned on and off by adjusting the po- 
tential depth V^i . If we first prepare each atom in the lat- 
tice into the superposition state (||) + |J.)) /\/2, and then 
lower V^i for a time T with A^ Z T = 7r/4 mod n/2, the fi- 
nal state is a cluster state with its dimension determined 
by the dimension of the lattice [11]. The rf-dimcnsional 
(d > 2) cluster states have important applications for 
implementation of scalable quantum computation with 
neutral atoms: after its preparation, one can implement 
universal quantum computation simply via a series of 
single-bit measurements only [11]. The use of such clus- 
ter states can significantly alleviate the stringent require- 
ments on separate addressing of the neighboring atoms 
in the proposed quantum computation schemes [17,18]. 
Although the present approach is somewhat slower that 
the cold collision scheme [17], it allows to take advantage 
of its simplicity and the reduced dephasing rate. 

As our second example, we explore the rich phase dia- 
gram of the Hamiltonian (2) in the presence of magnetic 
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FIG. 1. Phase diagrams of the Hamiltonians (2) for bosonic 
and fermionic atoms: (a) at zero magnetic field, (b) with a 
longitudinal field h z , and (c) (for bosons only) with a trans- 
verse field h x . Each phase is characterized by the following 
order parameter: I - z-Neel order; II - z-ferromagnetic order; 
III - ly-Neel order for fermionic atoms and a;?/-ferromagnetic 
order for bosonic atoms; IV and V - spin polarization in the 
direction of applied field, z and x respectively. 
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FIG. 2. 2a. The contours with the three potentials in the 
form of Eq. (5). The minima are at centers of the triangles 
when tpo — n/2. (2b). The illustration of the model Hamilto- 
nian (4). (2c). The schematic atomic level structure and the 
laser configuration to induce spin-dependent tunneling. 



fields. For simplicity, we assume a bipartite lattice and 
identical spin exchange constants for all links. Figure lb 
shows the mean-field phase diagram for bosonic particles 
in the presence of a longitudinal field h z . This diagram 
was obtained by comparing energies of the variational 
wavefunctions of two kinds: (i) the Neel state in the z 
direction {Si) — (— l) l e z ] (ii) canted phase with ferro- 
magnetic order in the xy plane and finite polarization 
in the z direction (o^) = e x cos9 + e z sin9. Here, 9 is a 
variational parameter, and e ZiX are unit vectors in the 
directions z, x. Transition between the z-Neel and the 
canted phases is a first order spin-flop transition [19] at 
h z = Z{\ 2 Z — A^) 1 / 2 [Z is the number of neighboring 
atoms of each lattice site), and transition between the 
xy-Neel phase and the z polarized phase is a second or- 
der transition of the XY type at h z = Z(X Z + Aj_). In the 
absence of transverse magnetic field one can use the ex- 
istence of two sub-lattices to change the sign of Aj_ using 
the transformation a^' v — > {—) l a^' v . Hence, fermionic 
atoms in the longitudinal magnetic field have the same 
phase diagram as shown in Figure lb, except that their 
canted phase has transverse Neel rather than transverse 
ferromagnetic order. Results of a similar mean- field anal- 
ysis of the Hamiltonian (2) for bosonic atoms with a 
transverse magnetic field h x are shown in Fig. lc. For 
fermionic atoms in a transverse field, there is one more 
phase with a Neel order along y direction. 

The third example involves the anisotropic spin model 
on a 2-D hexagonal lattice proposed recently by Kitaev 
[12]. In this model interactions between nearest neigh- 
bors are of the XX, YY or ZZ type, depending on the 
direction of the link 

H= £ A,<^, (4) 

u— x,y ,z\ 

where the symbol G D v denotes the neighboring 

atoms in the D v (y — x,y, z) direction (see Fig. 2b). 

To implement this model using ultracold atoms, we 
first raise the potential barriers along the vertical direc- 



tion Z in the 3-dimensional (3D) optical lattice so that 
the tunneling and the spin exchange interactions in Z di- 
rection are completely suppressed [9,1]. In this way, we 
get an effective 2D configuration with a set of indepen- 
dent identical 2D lattice in the X-Y plane. We then apply 
in the X-Y plane three trapping potentials (identical for 
all spin states) of the forms 

Vj (x, y) — Vb sin 2 [k\\ (x cos 9j + y sin 9j) + tpo] , (5) 

where j = 1,2,3, and Q\ = n/6, 62 = 7r/2, 6*3 = — tt/6. 
Each of the potentials is formed by two blue-detuned 
interfering traveling laser beams above the X-Y plane 
with an angle ^11 = 2 arcsin (1/V3) , so that the wave 
vector fcy projected onto the X-Y plane has the value 
fe|| = fcsin (<p\\/2) = k/V3. We choose the relative phase 
if{) = 7r/2 in Eq. (5) so that the maxima of the three 
potentials overlap. In this case, the atoms are trapped 
at the minima of the potentials, which form a hexagonal 
lattice as shown by the centers of the triangles in Fig. 
2a. We assume that there is one atom per each lattice 
site, and this atom interacts with the three neighbors 
in different directions through virtual tunneling with a 
potential barrier given by Vo/4. 

In such a hexagonal lattice, we wish to engineer 
anisotropic Heisenberg exchange for each tunneling direc- 
tion (denoted by D Xl D y , and D z , respectively). To this 
end, we apply three blue-detuned standing-wave laser 
beams in the X-Y plane along these tunneling directions: 

Vv„ (x,y) = V va sin 2 [k (xcose' v + ysin^)] , (6) 

where v = x, y, z, and Q' x = —n/3, 9' y = ir, 9' z — 7r/3. In 
general, we require that the potential depth V va depend 
on the atomic spin state as 

K CT = K + |+U+| + K-|-U-|, (v = x,y,z), (7) 

where 1+)^ (|— )„) is the eigenstate of the corresponding 
Pauli operator u v with the eigenvalue +1 (—1). 
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The spin-dependent potentials in the form of Eqs. (6,7) 
can be realized, for instance, with the specific atomic level 
configuration shown in Fig. 2c. Here, a =T,i denote 
two hyperfine levels of the atom with different energies. 
They are coupled to the common excited level |e) with 
a blue detuning A, respectively through the laser beams 
L\ and Li with frequencies matching the corresponding 
transitions. The quantization axis is chosen to be per- 
pendicular to the X-Y plane, and the phase-locked laser 
beams L\ and Li are both polarized along this direction. 
In the tunneling direction D z , we only apply the L\ laser 
beam, which induces the potential V za (x, y) with the de- 
sired form (7) of its depth V z<7 . In the tunneling direction 
D x or D y , we apply both lasers L\ and Li, but with dif- 
ferent relative phases, which realize the desired potential 
depth V xa or V yiJ of the form (7) in the corresponding 
direction. 

The potentials (6,7) do not have influence on the equi- 
librium positions of the atoms, but they change the po- 
tential barrier between the neighboring atoms in the D v 
direction from Vq /4 to Vl a = Vq/A + V ua . The param- 
eters V„+ and V"„_ in Eq. (7) can be tuned by varying 
the laser intensity of L\ and Li in the D u direction, and 
one can easily find their appropriate values so that in 
the D v direction, the atom can virtually tunnel with a 
rate t+ v only when it is in the eigenstate |+)„. Hence, 
it follows from Eqs. (2,3) that the effective Hamiltonian 
for our system is given by Eq.(4) with A„ s=a —t\ v / (2C7) 
for bosonic atoms with ~ f/f « ~ U. After com- 
pensating effective magnetic fields, we find exactly the 
model described by the Hamiltonian (4). 

The model (4) is exactly solvable due to the existence 
of many conserved operators, and it has been shown to 
possess very interesting properties [12]. In particular, it 
supports both abelian and non-abelian anyonic excita- 
tions, depending on the ratios between the three param- 
eters \ v . In the region where 2A„/ (A^ + X y + \ z ) < 1, 
(v = x, y, z), the excitation spectrum of the Hamiltonian 
(4) is gapless, but a gap opens when perturbation mag- 
netic fields are applied in the x, y, z directions, and the 
excitations in this case obey non-abelian fractional statis- 
tics. Out of this region, except at some trivial points 
with X x \ y \ z = 0, the Hamiltonian (4) has gapped exci- 
tations which satisfy abelian fractional statistics. Thus, 
the present implementation opens up an exciting pos- 
sibility to realize experimentally the exotic abelian and 
non-abelian anyons. 

Now we briefly discuss the techniques for probing the 
resulting states. To detect the quantum phase transi- 
tions in the XXZ model with magnetic fields or in Ki- 
taev's model, one can probe the excitation spectra via 
Bragg or Raman spectroscopy. In general, different quan- 
tum phases are characterized by specific dispersion re- 
lations (for instance, in Kitaev's model, one phase is 
gapped while the other is gapless). If the two probe 
light beams have the momentum and frequency differ- 



ences which match those of the dispersion relation in the 
corresponding phase, a resonant absorption of the probe 
light could be observed [21]. The direct observation of 
the fractional statistics in Kitaev's model can be based 
on atomic interferometry with a procedure similar to that 
described in Ref. [13]: one generates a pair of anyonic ex- 
citations with a spin-dependent laser focused on two lat- 
tice sites, rotates one anyon around the other, and then 
brings them together for fusion which gives different re- 
sults depending on the anyonic statistics. Other methods 
for detecting complex quantum states of atoms have also 
been developed recently [22]. 

In summary, we have described a general technique to 
engineer many-body spin Hamiltonians. 
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